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1. INTRODUCTION 



Since the early 1970s, many different techniques have been used to analyze the 
behavior of an antenna with a nonlinear load (References 1 through 11). While the 
characteristics of antennas with nonlinear loads have been analyzed in both the frequency 
and time domains, until now little attention has been given to possible chaotic effects that 
can result as a consequence of such nonlinearities. This oversight is not surprising, since 
only recently the field of chaos has moved from an essentially academic phenomenon to 
an intense applied research area in which chaotic behavior can be used to solve challenging 
technical problems. 

In the following discussion, an electrically small dipole antenna was loaded with the 
nonlinear circuit known as Chua's circuit (References 12 through 17). A particularly 
simple and very widely studied real nonlinear dynamical system, Chua's circuit has a 
number of significant advantages over other similar nonlinear circuits. First, it is the 
simplest circuit known (consisting of a linear inductance, a linear resistance, two linear 
capacitances, and only one nonlinear element) that still exhibits a rich variety of chaotic 
phenomena. Second, it is a physical system for which the presence of chaos has been 
proven analytically, simulated numerically, and demonstrated experimentally 
(References 14 and 17). Third, Chua's circuit is easily constructed using standard 
electronic components at low cost (Reference 18). It is an ideal example of a low-order 
physically simple system that can exhibit extremely complex, nonperiodic, bounded 
behavior. 

The reason for loading a dipole with Chua's circuit is to attempt to create an antenna 
that can easily switch from chaotic to nonchaotic behavior and vice versa. As is shown, 
this switching is easy to do in terms of the input voltage and current, as functions of time 
at the antenna input terminals, by simply changing certain values of the linear reactive and 
resistive elements (i.e., using variable resistors) that make up Chua's circuit. For an 
electrically small antenna (small in spatial extent such that the antenna length is only a 
small fraction of the wavelengths considered), it can be demonstrated numerically that 
when the antenna voltage function is temporally chaotic with a characteristic dynamical 
behavior the antenna will radiate temporally chaotic electromagnetic fields (referenced to 
some later time), exhibiting the same characteristic dynamical behavior. Similarly when 
the antenna voltage function is not chaotic, the associated fields (referenced to some later 
time) will be nonchaotic also. The ability to radiate a waveform that is predictable and 
periodic in time and then to suddenly switch to radiating a waveform that is complex and 
nonperiodic in time has obvious application to the area of secure communication. 
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In Section 2, the simplest equivalent circuit model for the electrically small dipole, 
modeled as a pure capacitance and loaded with a modification of Chua's circuit, known as 
Chua's oscillator, is shown and discussed. In Section 3, the relationship is derived 
between the voltages and currents resulting from the nonlinear load and the current that 
flows on the antenna element. Also an approximate vector potential function is 
introduced (based on the electrically small assumption) and the electromagnetic field 
components are determined. At this point both time series plots and delay coordinate 
plots of the electromagnetic field components are used to demonstrate when the 
electromagnetic fields are chaotic and when they are not. Section 4 provides discussion 
and conclusions. 



2. EQUIVALENT CIRCUIT OF A SMALL DIPOLE LOADED WITH CHUA'S 
OSCILLATOR: ANTENNA CAPACITANCE 



The original isolated Chua's oscillator, shown in Figure 1, is composed of two linear 
resistors, R = 1/G and Ro; two linear capacitances, C\ and C 2 ; a linear inductor, L; and a 
voltage-controlled nonlinear resistor, N R , called a Chua's diode (Reference 18). When R 
is set equal to zero, this circuit is known as Chua's circuit. Figure 2 shows the electric 
dipole loaded with Chua's oscillator. Imposing the restriction that the antenna is thin and 
electrically small and can be modeled as an induced voltage, V a (t), at the antenna terminals 
and an antenna capacitance, C a , (Reference 1), the antenna parameters and Chua's 
oscillator can be combined as in the equivalent circuit of Figure 3. In Figure 3, V a and C a 
are now in parallel with C 2 , V 2 and C b Vj. This equivalent circuit results in three coupled 
ordinary differential equations given by 



(la) 



C 2 + C L 



[h-(Va-Vi)G] 



(lb) 



•a 



(lc) 
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FIGURE 1. Chua's Oscillator. 
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FIGURE 2. Nonlinearly Loaded Dipole Antenna. 
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FIGURE 3. Equivalent Circuit of a Small Dipole Antenna Loaded With 
Chua's Oscillator. 
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Equation 1 describes a three-dimensional phase space (or state space) given by 



X(r) = {/ 3 (0,^(0,V1(0} 



(2) 



which is the state of the system at time t. A solution X(t) starting from some initial state, 
{I 3 (t ), V a (to), Vi(to)}, at some beginning time, t , is called a trajectory of Equation 1 . In 
Equation lc, f(V0 denotes the piecewise-continuous linear approximation of the driving- 
point characteristic of the Chua's diode. The function f^) is given by Reference 17, 



where E > 0, G a < 0, and G b < (Reference 17, Figure 3). 

Several explanatory remarks are in order. In Figure 3, we used the fact that V a = V 2 '. 
Also it is important to note that the antenna capacitance, C a , is usually on the order of 
10" 12 farads (F) for an electrically short thin dipole. This value is generally much less than 
typical values of C 2 (see Figure 3), which are usually on the order of 10" 7 F 
(Reference 19). Thus, in general C a « C 2 and C a « C } also. From further consideration 
of Equations 1 and 3, if the C 2 of Figure 3 is reduced by the exact amount, C a , then, when 
C a is added back again into Equation 1, this three-dimensional system is numerically 
identical to the isolated Chua's oscillator equations. An analysis of this set of equations 
is given in Reference 1 9. 

Once V a (t), the voltage at the antenna input terminals, is known from solution of 
Equation 1, the current on the antenna at the input terminals is given by 



/(viH 



'G b V x +(G b -G a )E ; V x <-E 
G a Vi ; M<|£| 

G b V { -{G b -G a )E ■ V X >E 



(3) 



I a (0 = C a V a (t) 



(4) 



or 




{ho-iVaV-mp} 



(5) 
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For the dynamical system of Figure 3, the three local Lyapunov exponents are given 

by 

h=~^ (6a) 
h = -7-777- (6b) 

L a + ( -2 



{G + G b ) 
(G + G a ) 



; outer regions 

(6c) 

; inner region 



It is obvious that X 3 is the only local Lyapunov exponent with the possibility of 
being positive (and thus producing chaotic behavior) (References 20 and 21) in an average 
sense. The sign of ^ 3 is, of course, dependent on the magnitude relationships among the 
linear conductances G, G a , and G b . Using Equation 1, and either Equation 4 or 5, the 
antenna current and voltage at the input terminals are known numerically as functions of 
time and are used in the next section to obtain the electromagnetic fields of the dipole 
antenna loaded with Chua's oscillator. 



3. VECTOR POTENTIAL FORMULATION 



Assuming that the small dipole of length 2h is oriented along the z-axis of a right- 
handed coordinate system, the associated vector potential can be written as 



h h 



Z',t- 



r-zz 



dz' 



4k 



-h 



\r-zz\ 



(7) 



At this point, if the antenna is small, one can assume that the current does not vary 
spatially along its length, z'. This is essentially an infinitesimal dipole assumption when 
the current is assumed to be constant in a spatial sense. (A triangular or sinusoidal spatial 
distribution could be used also if desired.) In this case, 



z'.t- 



\r-z'z^ 



\r-z'z\^ 



(8) 
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For an electrically small antenna, the first-order approximation 



r = r-z z 

can be used, and thus the right-hand side of Equation 8 becomes 



(9) 



j 



t-- 

\ c 



(10) 



Now that Equation 10 can be used to replace the current function in Equation 7 and 
is no longer dependent on z', the vector potential becomes simply 



cj 



(11) 



Now, I a (t) is already known from solving Equation 1 numerically and subsequently 
substituting the results into either Equation 4 or 5, and thus the time-retarded function in 

Equation 11 is obtained by substituting the new argument, t - into the known I a (t). 

c 

The new retarded time, t' = t — , is simply an indication that as the dipole radiation 

c 

moves outward in a spherical wave (at least, in the far-field) traveling at the speed of light, 
the radiation at time t' will not reach some fixed field point, r , until some later time, t. 

Using 



H = —V x A ; ^ = c 2 V xV xA 
HO dt 



(12) 



the electromagnetic field components in the time domain at some field point, r , and some 
time, t, can be written in the form 



2m- 



1 dV a (t') 1 1 d 2 V a (t') 
r dt' c dt' 2 



(13a) 



^^CJrcose 



EqKv 



r c dt' 



(13b) 



E e {T,t) = ^ hsin6 



£ 27ir 



idV a (t') i V a (t') { l d 2 V a {t') 
rc dt' r 2 c 2 dt' 2 



(13c) 
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Equation 13b and c was obtained by substituting Equation 4 into Equations 11 and 
12 and integrating indefinitely over the retarded time, t'. For now, the constants of 
integration have been set to zero for convenience. It is worth noting that in computing 
Equation 13, extensive use was made of the relationship 

%a(t') = l^a(Q = Idlgjt') 
dr c dt c df } 

Because the electromagnetic field components in Equation 13 can all be expressed in 
terms of V a (t') and its time derivatives, we speculate that if the antenna is sufficiently 
small for the equivalent circuit model of Figure 3 to be accurate and if the antenna input 
voltage is chaotic then the electromagnetic field components will be temporally chaotic 
also (at some later time). Similarly, if V a (t') is not chaotic, then the fields will not be 
chaotic at some later time. 



At this point it is necessary to demonstrate this translation of chaotic (or nonchaotic) 
behavior from the surface of the antenna, i.e., from the input terminals to the radiated 
electromagnetic field components. Usually time series plots of the electromagnetic field 
components (i.e., the field component amplitude plotted as a function of time) can be 
used. When a time series plot exhibits periodic behavior, the demonstration of periodicity 
is very compelling. However, chaotic states are not as easily shown to be chaotic simply 
from using time series plots. For possible chaotic states, it is more useful to employ a 
delay coordinate plot (Reference 22), which can reveal the characteristic shape of a 
chaotic attractor, if one exists. Thus, in Section 2 for X(t) given in Equation 2, it is 
possible to plot I 3 , V a , and Vi in a three-dimensional state space and to determine whether 
I 3 , V a , and V! exhibit chaotic behavior. Each variable could also have been plotted as a 
function of time in a time series, although this method would not have been as effective 
for demonstrating chaos. 

For the electromagnetic field components in Equation 13, these are known only as 
functions of time and can be plotted in time series also. But the delay coordinate method, 

for example, plots E e (rj) vs. Eg rj-- vs. Eq r,f . In this case, each value of 



cj 



c J 



some time series, i.e., Eg(r,t) vs. t, is plotted versus a time-delayed version of itself. If 
the delay coordinate plot closes on itself, i.e., if the orbits are closed orbits, then periodic 
behavior has been demonstrated. If the delay coordinate plot is very complex, not closing 
on itself and producing some type of strange attractor, then chaotic behavior has been 
demonstrated. Although delay coordinates are generally used to reconstruct an initially 
unknown state space from a limited set of measured data, they can also be used for 
functions such as Equation 13. 



9 



NAWCWD TP 8460 



To demonstrate numerically whether or not a chaotic (or periodic) voltage and 
current on the antenna will produce the same type of chaotic (or periodic) electromagnetic 
field behavior, three cases are considered: 

1 . A periodic antenna input terminal voltage, V a (t'), is used where the conductance, 

G = ~R 9 m Figure 3 is 530,0 microsiemens (nS). (All other values are given in the figure 
legends.) 

2. An antenna input terminal voltage, V a (t') 5 is a limit cycle of period-four, 
obtained when the conductance is increased to 539.0 |iS. 

3. A fully chaotic antenna input terminal voltage, V a (t') ? is a spiral Chua's chaotic 
attractor, obtained when the conductance is increased to 550.0 jllS. 

Thus, the route to chaos that has been chosen is the period-doubling route and the 
bifurcation parameter is G, the conductance (References 19 and 21). We demonstrate 
numerically, using both time series plots and delay coordinate plots, that the dynamics 
occurring at the surface of the antenna also occur at some later time in the radiated 
electromagnetic field. 



CASE 1: PERIODIC CASE 

The first case assumes that the linear circuit parameters of the equivalent circuit in 
Figure 3 have values such that V a (t') is periodic (Reference 19) (Figure 4a). Figure 4a is a 
three-dimensional phase space plot of I 3 , V a , and V b determined by solving Equation 1. 
This plot is a limit cycle of period- 1 and exhibits only one closed loop, immediately 
indicating the periodic nature of V a (t') in this case. Note that all circuit parameters are 
given in the figure captions. Figure 4b is the corresponding time series plot of the antenna 
voltage amplitude versus time and is oscillatory with a single period. 

Figure 5a is a three-dimensional delay coordinate plot of the magnetic field, H r at 
three separate times and separate fixed radii. Using Equation 13 and the extra parameters 

( r\ ( 2r\ 

shown in Figure 5a, HJt) vs.HJt — vs. t is plotted, where r is the radial 

Y \ cj Y \ c J 

distance to some field point away from the antenna and c is the speed of light. This delay 
coordinate plot of H^, while somewhat different in shape and orientation from the antenna 
voltage in Figure 4a, exhibits exactly the same dynamical behavior, i.e., a limit cycle of 
period-one. 
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(a) Three-dimensional phase space plot of antenna voltage versus Vi and I 3 . Circuit 
parameters are R = 12.5 Q,L= 18.0 mH, C, = 10.0 nF, C 2 = 100.0 nF, 
G a = -757.576 uS, G b = -409.090 uS, E = 1.0 V, C a = 10.0 pF, G = 530.0 uS. 



V a (t) 




(b) Time series plot (same circuit parameters as Figure 4a). 
FIGURE 4. Plot of Antenna Voltage Showing Periodic Behavior. 
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(a) Three-dimensional delay coordinate plot, H,,,. Circuit parameters are R = 12.5 Q, 
L = 18.0 mH, d = 10.0 nF, C 2 = 100.0 nF, G a = -757.576 uS, G b = -409.090 uS, 
E= 1.0 V, C a = 10.0 pF,G = 530.0 uS,r = 7.5x 10 3 m,h = 0.01 m,e = 7t/2, 

I = 100 A,c = 3x 10 8 m/s. 



H,(t) 




(b) Time series plot of magnetic field versus time (same parameters as Figure 5a). 
FIGURE 5. Plot of Magnetic Field Showing Periodic Behavior. 
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In Figure 5b, the periodic nature of the magnetic field as a function of time, is shown 
in a time series plot. Similarly, Figure 6a is a three-dimensional delay coordinate plot of 



K 



the electric field E e (since is assumed to be ^> the E r component is zero), at three 



separate times and separate fixed radii, i.e., E e (t) vs.Eq 



( A 
t-~ 

V c) 



2r\ 
VS. Eq\ t 



It too 



c J 



exhibits period-one limit cycle behavior just as shown in Figure 4a. Figure 6b is the 
corresponding time series plot for E as a function of time and is periodic. 



E e t- 



0.002 



-0.002 



0.002 



EJt- 



2r 



-0.D02 




(a) Three-dimensional delay coordinate plot, E e , (same parameters as Figure 5a). 



FIGURE 6. Plot of Electric Field Showing Periodic Behavior. 
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E e (t) 




(b) Time series plot of electric field versus time (same parameters as Figure 5a). 

FIGURE 6. (Contd.) 

CASE 2: PERIOD-4 CASE 

The second case assumes that the parameters of the equivalent circuit in Figure 3 are 
the same as in Case 1, with the exception that now G = 539 |uS. Increasing the value of 
the conductance causes the limit cycle of Case 1 to now bifurcate so that it closes on itself 
after four loops of the trajectory. This is a result of period-doubling, a well known route 
to chaos (References 19, 21, and 22). This intermediate case is shown in Figure 7a, where 
the phase space dynamics involving the antenna voltage, V a , exhibit four loops before 
closing on itself, indicating period-4 limit cycle behavior. The four distinct periods are 
also shown in the time series plot of V a versus t in Figure 7b. 

Using the same parameters as those used in Figure 7a, Figure 8a is a three- 
dimensional delay coordinate plot of the magnetic field, which also shows four separate 
loops indicating period-4 behavior. Although amplitude, orientation, and loop shape are 
different from those in Figure 7a, the period-4 dynamics are undeniable. Figure 8b is the 
corresponding time series plot with four distinct periods. Similarly in Figure 9a and b, the 
delay plot and the time series plot for E e exhibit the same period-4 behavior. 
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(a) Three-dimensional phase space plot of antenna voltage versus V! and I 3 . 
Circuit parameters are Rq = 12.5 Q, L = 18.0 mH, C x = 10.0 nF, C 2 = 100.0 nF, 
G a = -757.576 uS, G b = -409.090 uS, E = 1.0 V, C a = 10.0 pF, G = 539.0 uS. 



v a (t) 
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Time series plot of antenna voltage versus time (same circuit parameters as Figure 
FIGURE 7. Plot of Antenna Voltage Showing Period-4 Behavior. 
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(a) Three-dimensional delay coordinate plot of magnetic field, H,,,. Circuit parameters 
Ro = 12.5 Q, L = 18.0 mH, C, = 10.0 nF, C 2 = 100.0 nF, G a = -757.576 uS, 
G b = -409.090 uS, E = 1.0 V, C a = 10.0 pF, G = 539.0 uS, r = 7.5 x 10 3 m, 
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(b) Time series plot of magnetic field versus time (same parameters as Figure 8a). 
FIGURE 8. Plot of Magnetic Field Showing Period-4 Behavior. 
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(a) Three-dimensional delay coordinate plot of electric field, E, 
(same parameters as Figure 8a). 



E e (t) 




Time series plot of electric field versus time (same parameters as Figi 
FIGURE 9. Plot of Electric Field Showing Period-4 Behavior. 
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CASE 3: CHAOTIC CASE 

The third case assumes that the circuit parameters of the equivalent circuit of 
Figure 3 are the same as for the two previous cases, except that now the conductance has 
been increased to 550 \iS (Reference 19). This particular value of G causes the antenna 
voltage, V a (t'), to be chaotic, specifically resulting in a spiral Chua's chaotic attractor 
(Reference 19). It is shown in Figure 10a and is an extremely distinctive chaotic attractor. 
Figure 10b is the associated time series plot of the antenna voltage amplitude as a function 
of time, which immediately shows complex, bounded, nonperiodic behavior. 

Figure 1 la is the three-dimensional delay plot of the magnetic field, H^. While again 
showing different amplitude, shape, and orientation when compared to Figure 10a, it 
exhibits the same spiral chaotic attractor form, including the central loop that overlaps the 
remainder of the spiral, with the remaining overlapping loops occurring at the top of the 
figure rather than at the bottom. A comparison of Figure 10b and Figure lib shows 
similar behavior in the time series plots. Figure 12a, the delay plot of E e , also gives a very 
faithful representation of the spiral Chua's chaotic attractor, again showing the same 
overlapping loops as in Figure 10a in somewhat different orientations but distinctively 
seen. The time series plot of E e in Figure 12b is almost identical in form to the time series 
antenna voltage plot in Figure 10b. 
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(a) Three-dimensional phase space plot of antenna voltage versus Vj and I 3 . 
Circuit parameters are Rq = 12.5 Q, L = 18.0 mH, Q = 10.0 nF, 
C 2 = 100.0 nF, G a = -757.576 uS, G b = -409.090 uS, E = 1.0 V, C a = 10.0 pF, 

G = 550.0 uS. 



v a (t) 




(b) Time series plot of antenna voltage versus time showing complex nonperiodic 
behavior (same circuit parameters as Figure 10a). 

FIGURE 10. Plot of Antenna Voltage Showing Chaotic Behavior. 
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(a) Three-dimensional delay coordinate plot of magnetic field, H,,,. Circuit parameters are 
Ro = 12.5 Q, L = 18.0 mH, C, = 10.0 nF, C 2 = 100.0 nF, G a = -757.576 uS, 
G b = -409.090 uS, E = 1.0 V, C a = 10.0 pF, G = 550.0 uS, 
r = 7.5x 10 3 m,h = .01 m,6 = 7t/2, I a = 100 A,c = 3x 10 8 m/s. 



H $ (t) 




(b) Time series plot of magnetic field versus time showing complex nonperiodic behavior 

(same parameters as Figure 1 1 a). 

FIGURE 11. Plot of Magnetic Field Showing Chaotic Behavior. 
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(a) Three-dimensional delay coordinate plot of electric field, E, 
(same parameters as Figure 1 la). 



E e (0 




Time series plot of electric field versus time showing complex, nonperiodic behavior 

(same parameters as Figure 11a). 

FIGURE 12. Plot of Electric Field Showing Chaotic Behavior. 
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4. CONCLUSION 



The electrically small dipole antenna loaded with Chua's oscillator as the nonlinear 
load was analyzed. By using certain numerical values for the linear reactive and resistive 
elements of this circuit, various types of periodic, period-doubled, and chaotic behavior 
that occur for the voltage at the antenna input terminals can also be shown to occur for 
the associated radiated electromagnetic field components. Using three cases, i.e., periodic, 
period-4, and chaotic behavior, we show that the voltage characteristic dynamics 
determined at the input terminals are the same dynamics of the radiated electromagnetic 
field components. 

Thus, for an electrically small dipole, we demonstrated that when the antenna voltage 
function is temporally chaotic, the antenna will radiate temporally chaotic electromagnetic 
fields with the same dynamical behavior referenced to some later time. Similarly, when 
the antenna voltage is periodic, the radiated fields will be periodically referenced to some 
later time also. Thus, switching between a periodic and a chaotic temporal 
electromagnetic field behavior can be accomplished by simply increasing the conductance 
of a single resistor of the antenna/Chua's oscillator equivalent circuit. 
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